Universal scaling in a trapped Fermi super-fluid in the BCS-unitarity crossover 
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Using numerical simulation based on a density-functional equation for a trapped Fermi super-fluid 
valid along the BCS-unitarity crossover, we establish robust scaling over many orders of magnitude 
in the observables of the system as a function of fermion number. This scaling allows to predict 
the static properties of the system, such as energy, chemical potential, etc., for a large number of 
fermions, over the crossover, from the knowledge of those for a small number (~ 4 — 10) of fermions. 
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I. INTRODUCTION: 

The Pauli principle among identical fermions leads to 
an effective repulsion which avoids the collapse of the 
hadronic universe. A many-fermion system dominated 
by Pauli repulsion exhibits universal behavior with the 
physical attraction playing a minor role. Such universal 
behavior manifests in correlation and scaling in the ob- 
servables of few- and many-fermion systems [l| . This uni- 
versality is prominent in the limit of zero Fermi- Fermi in- 
teraction in the Bardeen-Cooper-Schrieffer (BCS) theory 
of super-fluid (SF) fermions leading to universal prop- 
erties of low-temperature superconductors 0, of cold 
neutron matter and neutron star and of a trapped 
Fermi SF (3] at zero temperature. This universality also 
manifests at unitarity 0, 0], when the atomic scatter- 
ing length a tends to infinity (a —f oo). Both at the 
BCS limit (a — > —0) and unitarity, the scattering length 
ceases to be a scale of length, and the only length scale 
of the problem in these limits is kp^ where hkp is the 
Fermi momentum. Consequently, all energies are uni- 
versal functions of Fermi momentum kkp (or of Fermi 
energy Ep — h^kp/2ni), where m is the atomic mass, 
and the system acquires universal behavior with robust 
scaling involving the stationary observables, such as en- 
ergy, chemical potential, etc., as a function of number of 
fermions. 

How the above scaling appears in a fully paired N- 
fermion SF with equal number of spin up and down atoms 
confined in a spherically-symmetric trap U{r) can be seen 
in the local-density approximation (LDA) [4|, i2, iS] , where 
the system obeys 



C/(r) + ^[3^^n(r)m 



(1) 



where ^ = 1 in the BCS limit and = 0.44 (obtained Q 
from Monte Carlo calculations on a uniform Fermi SF) 
at unitarity, n(r) is the SF density, fxo is the chemical 
potential, and U(r) ~ muj'^r'^ /2 is the harmonic trap with 
Lo the trap frequency and normalization j n{v)dv ~ N. 
The second term on the left-hand side of Eq. ([T]) is the 



chemical potential of a uniform SF Fermi gas. The total 
chemical potential /i = I-J^qN and energy E of the traped 
system are given by 



11 = N J dr 
E = N dr 



2 m 

C/(r) + ||^(37rMn 
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nir), (2) 
n(r). (3) 
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It is useful to work in terms of dimensionlcss variables 
in harmonic oscillator units, obtained by setting ft = 
m = Lo = I, when lengths are expressed in units of 
I = ^ h/{mijj) and energy in ftw. The normalization con- 
dition for the density n leads to \/n{r) = [2(3iV)^/'^v^~ 
r2]^V2/(37r2) for 2{2,NY'^^ > and otherwise. Us- 
ing this in Eqs. © and © we get ^ = (3A^)'*/3^/3, 
E = {3N)'^^^^/^/4:. An extended LDA (ELBA) was sug- 
gested for this energy to include corrections for small 
N: E = {3N)^/^y/^{l+a/N^/^)/A, a = 0.5 {a chosen 
to provide agreement with Monte Carlo data d, 0] for 
trapped Fermi SF). 

The energy E = {SN^/^/A in the BCS Hmit (^ = 1) 
also follows from a consideration of the A'^-fermion sys- 
tem in a spherical harmonic trap at zero temperature. 
The degeneracy of the ith state of harmonic oscillator 
of energy haj{i + 3/2) is {i + + 2)/2. If we consider 
many paired states up to i = /(/ >> 1) completely full, 
then Ep = / (in units of fko) and the total number of 
fermions, considering two (spin up and down) fermions 

in each state, is TV = ELo(* + !)(« + 2) = (/ -h !)(/ + 
2)(/ + 3)/3 ~ 1^/3. The total energy in oscillator units is 
E = ELo*(* + 1)(« + 2) = /(/+l)(/+2)(/+3)/4 « /V4. 
Consequently, Ep = {3N)^/^,kp = V2{3N)^/^, and 
E = {3N)'^/^/4:, consistent with the LDA above. The 
energy at unitarity E = {3N)'^^^y^/4: differs from this 
energy at the BCS limit by a numerical factor to take 
into account the atomic attraction. The chemical poten- 
tial, with dimension of energy, differs from energy by only 
a numerical factor consistent with the universal behavior. 

The E and ^ vs. N scalings valid in the BCS and uni- 
tarity limits are useful in predicting E and fi for a large 
N from a knowledge of those for a small N. In this pa- 
per we address two interesting questions related to these 
scalings. (a) Are the scalings dependent on the LDA? (b) 
Can these scalings valid at the BCS and unitarity limits 
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be extended along the BCS-unitarity crossover (loj? 

The LDA, lacking a gradient term, does not include 
properly the variation of the SF density near the surface. 
This can be remedied by, following a suggestion of von 
Weizsacker , including a proper gradient term in the 
LDA and thus transforming it into a density-funcional 
(DF) equation Corrections to the gradient term have also 
been considered Q . The coefficient of the gradient term 
is fixed by requiring the DF equation to be Galilei in- 
variant and also equivalent to the hydrodynamical flow 
equations [T^ . The crossover from the weak-coupling 
BCS limit to unitarity has been an active area of 
research 0, [^, [l^ after the experimental realization 
of this crossover in a trapped dilute Fermi SF by ma- 
nipulating an external background magnetic field near a 
Fcshbach resonance. Also, by modifying the bulk chem- 
ical potential term in the LDA to one appropriate for 
the BCS-unitarity crossover, the above scalings can be 
extended to the full crossover. 
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FIG. 1: (online color at www.lphys.org) Energy and chemical 
potential of a trapped Fermi SF in oscillator units as a func- 
tion of atom number at unitarity obtained from a solution of 
Eq. © and from FNMC ^ and GFMC [l] calculations as 
well as LDA and ELDA calculations. 



II. DENSITY-FUNCTIONAL FORMULATION: 

To study the universality along the crossover, we use a 
Galilei-invariant DF formulation of the trapped Fermi SF 
[T^ . [isl . flGj , equivalent to a hydrodynamical model with 
the correct phase- velocity relation 3 v = hV9/{2m), 
where v is the SF velocity, and the phase of the order 
parameter at position r: 



8m 
/i(n,a) 



+ U + ^i{n,a) 

ft2 



2m 



(37r2n)2/3g(„i/3„)^ 



Mo 



g{x) = 1 + 



fi = N I dr 
E = N I dr 



{1 - f3ix + I32X^) 
8m 



|V\/n|^ + Un + fi{n,a)n 



„ \S/y/^f + Un+ n{n\a)dn' 
8m Jo 



(4) 

(5) 
(6) 
(7) 
C8) 



where XItX^tPi a-nd P2 are parameters. If we choose 
XI = 47r/(37r2)2/3,^2 = 300, /3i = 40, /32 = X2/(l - 
the DF equation ^ (i) produces energies that agree 
with the fixed-node Monte Carlo (FNMC) [i| and Green- 
function Monte Carlo (GFMC) Q results for the energies 
of a trapped Fermi SF at unitarity (for N < 30), and 
over the crossover [HI, [13 for A^ = 4 and 8, (ii) provides 
a smooth interpolation between the energies of a SF at 
the BCS and unitarity limits [l^, [T3l. (hi) Moreover, 
the bulk chemical potential u(n, a) of Eq. ^ is exact 
at unitarity and satisfies [T^ . [Tq the weak-coupling BCS 
limit lim^^o3(2;) I + inx/iSiT^)^/^ [Hi. 

The gradient term in Eq. (|4]) corresponds [l2| to 
a quantum pressure in the equivalent hydrodynamical 
equations and provides a correction to the LDA [31 . LDA 
is a good approximation for a large N, when the bulk 



chemical potential ^{n,a) — a positive term responsible 
for Pauli repulsion in the system even for attractive (neg- 
ative) a — is very large. The gradient term is consistent 
with the hydrodynamic flow of paired fermions of mass 
2m [l,[l3l. 



III. NUMERICAL RESULTS OF SCALING 

We solve the DF equation ^ by the split-step 
Crank-Nicolson method by transforming it into a time- 
dependent equation and using the FORTRAN programs 
provided in Ref. [T^. The space and time steps used 
in discretization of the equation were 0.05 and 0.001, re- 
spectively. 

The calculated energy and chemical potential of the 
trapped Fermi SF as a function of the number of Fermi 
atoms are exhibited in Fig. [TJ We also compare the 
results for energy with those obtained by the FNMC 
and GFMC ^ calculations as wcU as LDA and ELDA 
calculations. It is found that the DF and ELDA results 
for energy are in better agreement with the Monte Carlo 
simulations than the LDA results. 

An analysis of the DF results for energy and chemical 
potential of Fig. [T] reveals that these functions maintain 
the following quantities 

S{N) = [E{N)/N^/^ - 0.37l/A^2/^ (9) 
k{N) = [^{N)/N^/^ - 0.27l/A^2/^ (10) 

fixed at constant values independent of N. This is a 
consequence of scalings E{N)/N'^/'^ and ^{N)/N^^^ vs. 
N^^^ as shown in Fig. [T] If the functions S{N) and k{N) 
were really universal, then they should maintain approx- 
imate constant values along the BCS-unitarity crossover 
independent of N, although weakly dependent on the 
atomic scattering length a. To demonstrate it, we plot 
in Fig. ^E{N)/N^^^ and fi{N)/N^^^ vs. A^^/s 
tained from the DF equation ^ along the BCS-unitarity 
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FIG. 2: (online color at www.lpliys.org) Same as in Fig. [T] 
but with results for different scattering lengths a: a = — 
(BCS limit), a = —0.5 (in BCS-unitarity crosover), a = — oo 
(unitarity). (The LDA results are not shown.) 
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FIG. 3: (online color at www.lphys.org) (a) 5{N, kpa) — 
[E{N,kFa)/N'-2/3) - 0.37] /A*^' 2/3) vs. [kpal for differ- 
ent N from the solution of DF equation (|4| (lines) and 
from FNMC calculations ^ (points), (b) K{N,kFa) = 
[fi{N,kFa)/N^2/3) - 0.27]/Af(2/3) vs. \kFa\ for different N 
from the solution of DF equation ([4} (lines). 

crossover for three values of a representing the BCS limit, 
unitarity and one in the crossover. Figure [2] clearly illus- 
trates the robust scaling. 

To take the full advantage of the above scalings 
at unitarity, as shown in Figs. [T] and [2l along 
the full BCS-unitarity crossover we plot (5(-/V, fc^a) = 
[E{N,kFa)/N^/^ - 0.37]/iV2/3 and K{N,kFa) = 
[n{N,kFa)/N'^^^ - 0.27]/iV2/3 vs. Ikpa] in Figs. [3] (a) 
and (b), respectively, where kp = \/2(3Af)^/^. It should 
be noted that kpa is a model-independent dimension- 



less measure of atomic interaction and should hence be 
used in the study of universal scaling. If the Fermi SF 
were really dominated by the universal scaling, a plot of 
S(N, kpa) and k{N, kpo) vs. \kFa\ should lead to univer- 
sal curves indepndent of N and this is indeed so as can 
be seen from Figs. [31 

TABLE I: Dimensionless energies E/lhLuiZNf'^] of a trapped 
Fermi SF along the crossover. The last two columns, {4} 
and {10}, give the predicted energies for N — 10^ employing 
scaling (lllfl using the energies for A' = 4 and 10, respectively. 



a\N = 


10 


10^ 


10^ 


10^ {4} 


10^{10} 


-0.001 


0.2656 


0.2539 


0.2499 


0.2428 


0.2470 


-0.01 


0.2652 


0.2531 


0.2459 


0.2393 


0.2434 


-0.1 


0.2550 


0.2364 


0.2046 


0.2027 


0.2049 


-1 


0.2079 


0.1869 


0.1714 


0.1731 


0.1740 


-10 


0.1897 


0.1731 


0.1665 


0.1686 


0.1693 


-100 


0.1875 


0.1714 


0.1659 


0.1682 


0.1688 



The universal nature of the energy curves as illustrated 
in Fig. [3] (a) allows us to predict energy of the N fermion 
system from a knowledge of the energy of N' fermions 
from the universal relation 

E{N,a)/N^/^ - 0.37 _ E{N' , a)/N'^/^ - 0.37 

Ar2/3 ~ Ar'2/3 ' ' 

We have done so and our results for energy as calculated 
by solving the DF equation are exhibited in Table |T] [3 , 
where we display energy along the crossover for different 
N and a. We also calculated the energies for N = 10'^ 
atoms using the energies for A^' = 4 and 10. The pre- 
dictions so obtained listed in Table 1 , compare well with 
the calculated results. This demonstrates the usefulness 
of the scaling ([TT|) for energy. 

IV. SUMMARY AND DISCUSSION 

From a numerical study of the static properties of a 
trapped Fermi SF using a Galilei-invariant DF formula- 
tion [l^l, equivalent to a generalized hydrodynamic for- 
mulation with the correct phase- velocity relation [J, we 
establish that, because of the dominance of the Pauli re- 
pulsion, the trapped Fermi SF has a universal behavior 
not only in the BCS and unitarity limits but also along 
the BCS-unitarity crossover. This allows a prediction of 
the static properties of a large Fermi SF in the crossover 
region from a knowledge of the same of a small system 
through a scaling relation, cf. Eq. ([TTjl . The predicted 
energies of a Fermi SF with 10^ atoms from a knowledge 
of the same with 4 (10) atoms is found to have an er- 
ror of less than 3%. Although we used a DF equation 
in our study, the observed scaling(s) should be indepen- 
dent of the use of the DF equation and also of the LDA. 
The scaling(s) is a consequence of the dominance of the 
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Pauli principle leading to a repulsive interaction, where 
the actual physical attraction plays a minor passive role. 

In a Bose- Einstein condensate (BEC), there is no Pauli 
repulsion, and the physical interaction plays a major ac- 
tive role. Consequently, in a BEC, the stationary ob- 
servables are very sensitive to the atomic interaction and 
no scaling should exist in a stationary observable as a 
function of N from weak-coupling to unitarity crossover. 



Nevertheless, there have been systematic studies of static 
properties of a BEC as a function of N . Fermi-Bose 
mapping; fermionization, and other considerations of a 
strongly-interacting one-dimensional Bose gas are taken 
up in Ref. [2l|. 

The work was partially supported by the FAPESP and 
CNPq of Brazil. 
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